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Abstract 
This paper proposes two modified generalized Gauss-Seidel iteration techniques to determine the Absolute Value 
Equations (AVEs). Convergence of the new techniques is established under some appropriate conditions lastly; several 


numerical examples verify the significance of the techniques. 
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1 | Introduction 


COD iccnsce 


i In this article, we explore the AVE of the form; 
Computational 


oa pape ae Ax — |x| = b. (1) 
Numerical Dimensions. 

This article is an open 

access article distributed | Here A € R’",b,x € R” and |x| represents the vector with the absolute value of each element of x. 
under the terms and 


conditions of the Creative] The AE (7) has attracted significant interest in the optimization domain for more than twenty 


Commons Attribution years. The most important reason is that they are closely related to Linear Complementarity Problems 
(ope leave (LCP) [1], [2] and horizontal LCP [3], which include a vast range of mathematical programming 
(http://creativecommons. : SS 

a ave a ange of applicatio -[10]. 
croflishasalby/LO) problems and have a wide range of applications [4]-[10] 


In recent years, many efforts have been made to develop techniques for finding numerical solutions 
to the AVE (7). For example, Salkuyeh [11] proposed the Picard-HSS approach to compute Eg. (7). 
The SOR-like system is discussed in [12]. Additionally, Chen et al. [13] explored various convergence 
effects for a SOR-like strategy using different parameters. In [15], the authors presented the shift 
splitting technique to obtain Eg. (7) and established the convergence properties of the offered 
approach. In the study of Moosaei et al. [16], two approaches are offered to determine Eq (7), the 
Homotopy perturbation approach and the Newton approach with the Armijo step. 


Corresponding Author: rashidali0887@gmail.com 
& https://doi.org/10.22105/cand.2022.157605 


The modification of the generalized gauss-seidel iteration techniques for absolute value equations 


Pakharzadeh and Shams [17] offered a Mixed-Type Splitting (MTS) strategy for the calculation of Eg. (7) 
and demonstrated its convergence conditions. In addition, there are several other numerical techniques in 
the literature to determine Eg. (7) [18]-[24]. 


The present study extends the Generalized Gauss-Seidel (GGS) method [14] and gives Modified GGS 
(MGGS) techniques, which can be improved by adding three extra parameters. In addition, we examine 
the convergence of the novel techniques under appropriate situations. 


An overview of this paper is provided below. In Section 2, we discuss the auxiliary results. Section 3 
examines the proposed techniques and their convergence. In Sections 4 and Sections 5, we present 
numerical results as well as conclusions. 


2 | Preliminaries 


This section provides a helpful lemma and a _ few _ notations for later analysis. 
Throughout this paper, |A| = ({a;;)) and ||Al|.. represent the absolute value as well as the infinity norm of 
A = (a;) € R™", respectively. In addition, the symbol ‘I signifies the identity matrix. 


Lemma 1 ([31]). Let x,z € R". Then |||4|- |xlll. < [l¥ - lle: 
3 | Proposed Techniques 


Here, we briefly describe the MGGS techniques. We have divided this section into two sections. The first 
section consists of the MGGS I and its convergence, and the second section concerns the MGGS II and 
its convergence to solve the AVEs. 


3.1 | MGGS I for AVE 


Let we split the A matrix as 


1 
A=— M-N), (2) 


With 
M=Q+Ap-AA,, N=Q+ 1-B)Ap+ B-AAL+ PAu. 


Where B and A are real parameters with B # 0 and QO = @J where0 < 8 <1. Also, A;, Ap and Aj, are strictly 
lower-triangular, diagonal and strictly upper-triangular parts of A. Based on Eg. (2), we can write Eq. (7) as 
follows: 


1 
F M -N)x — |x| =b, 
(3) 
Mx — B|x| = Nx + Bb. 
By using the iterative scheme, we define MGGS | as follows: 
Mxkt1 — Back] =Nxk + Bb, 


(4) 
x1 = BM xk] + M-1(Nxk + Bb), k = 01,2,...). 


Remark. If 6 = A =1 and Q = 0, then the MGGS I tranform to the GGS approach [14]. 


The subsequent theorem specifies the convergence of MGGS I. 


Theorem 1. ConsiderA = M-N, which represents a splitting of A with entries in the diagonal of A 
greater than one. If 


Mom], <1 pI" I... 6) 


Then the {x} sequence obtained from the MGGS I converges to the unique solution x* of Eg. (7). 


Proof. The proof of |||} <1 is easy and is skipped here (see Theorem 3 in [14]). Let us consider x* 


and h* to be two distinct solutions of Eg. (1) for the purpose of determining the uniqueness of the 
solution. By applying Eg. (4), we determine 


x* = BM" |x*| + M1} Mx* + Bb). (6) 
h* = BM" |h*| + M+ Mh* + pb). (7) 
Based on Eg. (6) and Eg. (7), we can get 
x* —h* = BM |x*| -|h*|) + MIN x* -h*). 
According to Lemma 1 and Eg. (5), we have 
[e* =F" los S BMI lie ie lle MOEN] Ibe Sailer 
[le — Aloo < PIMA] le — FA lloo + (= BI )Ib* — 1 lho, 
lt = leg BM | lle = 1 les (= BMI )lbe* =H loos 
(1 - p||M- |) o* - Ihe < (1 pI) lhe - 1 lho, 
|e =" lly < [le = Laas 
The result is a contradiction. Consequently, x* = h*. 
To determine convergence, we assume that x* is a unique solution of Eg. (1). Accordingly, from Eg. (6) 
and 
xh = BM | 4 (Nx + pb). 
We deduce 
kt] _ yk = pM (|x| _ *1) + MN(x! = a”), 
Taking the infinity norm, we get; 
[et Ts BMT et 2 + MNT bE I 
[et — |] — BMT ee 21, < ION | Ibe 1 
(1 pila", — 21], < MINI Ib - 21, 


vay NL 
coo 1 


k+1 _ a 
[h = BIIM ho 


Leo 
This inequality implies that convergence of the approach can be achieved if Eg. (3) is fulfilled. 
3.2 | MGGS II for AVE 


This unit offers the MGGS II. Using Eg. (1) and Eg. (2), the MGGS II is expressed as follows: 


i 
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Mrk+1 — | ec | = Nyk+1 + bb, 


x1 = BM c+] 4. M1 (Nak + Bb), kaU12 a. 
To demonstrate the convergenceiof MGGS I, we apply the following theorem. 
Theorem 2. Consider A = M — N, which represents a splitting of A with entries in the diagonal terms of 


A greater than 1. Then the {xk} sequence obtained from the MGGS IJ converges to the unique solution 
ot dg. (7). 


Proof. From Theorem 1, the uniqueness follows directly. In order to verify convergence, we can consider 
the following: 


xkl x = BMo1 x1) 4 M-1(Nxk+1 + pb) — (BM~!|x*| + M-! Nx* + Bb), 


M(x‘+1 _ ie") a B([x**}| - x*]) + N(x _ x), 

(8) 
Mxkt1 = Nxk+1 = Blxk+1| = Mx* —-Nx* = BIx*|, 
(M — N)xk+? — plxk*1] = (M-N)x* - Bix*l. 


Since 
1 
ae M-N) > fA=M-N. 


As a result, Eg. (8) can be expressed as 
BAxk+1 = B]x**+1| = BAx* _ BIx*|, 


Axktl - peor = Ax* — Ix*|, 


Axktl bel =b. 


Thus, x**! solves the AVE (1) system. 
4 | Numericaliexperiments 


This part gives four examples in order to provide an overview of the performance of the newly formulated 
techniques from three points of view: 


I. The iteration steps (demonstrated as ‘Itr’). 
II. The CPU time in seconds (denoted as “Time’). 
Ill. The residual error (represented as ‘RSV’). 


The term ‘RSV’ is described as follows: 


_[ax-bl-bl, 


RSV: <10~°. 
I|bl|2 


Example 1. Let 


E Rm, 
-1 
-1 4 
Compute b = Ax* — |x*| € R", where x* = (-1)*. Here the starting iterate is zero vector, and we compare 
both the proposed techniques with the SOR-like approach outlined in [12] (shown as SL) and the special 
shift splitting approach described in [15] (shown as SS). Tab/e 7 describes the results. 


Table 1. A comparison of example 1 using@ = 0.1, B = 0.98 andA = 0.8. 


Methods n 1000 2000 3000 4000 
SL Ttr 18 18 18 18 
Time 3.0156 13.1249 33.9104 65.1345 
RSV 6.12e — 07 6.13e — 07 6.12e — 07 6.14e — 07 
SS Iter 14 14 14 
14 
Time 2.8128 9.0954 17.3028 29.1644 
RSV 8.91le-—07 8.92e-07 8.93e — 07 8.93e — 07 
MGGS I Ter 9 9 9 9 
Time 2.0823 5.3422 12.3724 21.2038 
RSV 9.86e — 07 6.97e — 07 5.69e — 07 4.92e — 07 
MGGS I Itr 7 7 7 rg 
Tim 0.9208 1.8422 2.3677 6.2190 
RS 1.82e — 07 1.83e — 07 1.83e — 07 1.83e — 07 


Example 2. Let A =1+M € R™" and b = Ax* — |x*| € R" with 


> =0,5! 
-151 § ee : 
M = ee ening = (=) (R=12,.9n) €R™ 
me) 
-15I §& 
Where S = tridiag -1.5,4,—0.5) € R™1,I € R™4 and n = q’. In example 2, the same starting guess and 
stopping criteria as in [17] are used. We compare the suggested techniques with the AOR procedure [25] 


and the MTS scheme [17]. We illustrate our results in Table 2: 


Table 2. A comparison of example 2 using@ = 0.2, = A = 0.9. 


Methods n 100 400 900 
1600 
Ter 97 190 336 706 
AOR Time 0.4721 2.8203 3.2174 6.3887 
RSV 9.80 — 07 9.6le — 07 9.73e — 07 9.84 — 07 
Ter 88 157 250 386 
MTS Time 0.4041 1.7953 3.0219 5.7626 
RSV 8.91e - 07 9.65e — 07 9.18e — 07 9.56e — 07 
Ter 46 70 91 
110 MGGSI Time 0.2108 0.3844 0.9534 1.8747 
RSV 8.83e — 07 8.94e — 07 817e-07 8.57e-07 
Ttr 23 34 45 55 
MGGS II Time 0.1722 0.2344 0.7603 1.4200 
RSV 5,08e — 07 9.96e — 07 7.8le-07 _7.5le-07 
Example 3. Let 
8 -l 
18 2. 
M= 6B Vemma. 
-1 
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And b = Ax* — |x*|,€ R", where x* = (-1)*. Based on the same initial guess and ending criteria as explained 
in [26]. We compare our newly developed techniques with those presented in [26] (signified as SI with w = 
1.0455) as well as the SOR-like strategy offered in [12] (considered as SOR). 


Table 3. A comparison of example 3 using@ = 0.2, B = 0.95 and A = 0.85. 


Methods n 1000 2000 3000 4000 
5000 

SI Itr 13 13 14 14 14 
Time 3.9928 8.8680 24.4031 51.3946 73.3394 
RSV 6.04e — 07 8.54e — 07 2.33e — 07 2.69e —07 3.01le — 07 
SOR Itr 12 13 13 al; 13 
Time 1.5136 3.3817 6.1262 7.1715 9.5262 
RSV 9.45e — 08 2.69e — 08 3.29e — 08 3.80e-08  4.25e —07 

MGGS I Itr 11 11 11 11 11 
Time 1.2718 2.3422 3.0127 4.4824 6.2511 
RSV 9.45e — 08 2.69e — 08 3.80e — 08 3.80e -07 3.80e — 07 

MGGs I Itr 5 6 6 6 6 

Time 0.45532 0.5203 0.9376 1.2788 3.7826 
RSV_ 9.29e — 07 1.22e — 08 1.47e — 08 1.68e — 08 3.8le — 07 


Based on the results of Table 1 to Table 3, all tested approaches can compute the solution to Eq. (7). 
Nevertheless, we notice that both the ‘Itr’ as well as the ‘time’ of the suggested techniques are better than 
the other known procedures. To conclude, we can state that our presented techniques are practicable and 
effective for AVEs. 


5 | Conclusion 


We have shown two MGGS techniques for determining AVE (7) and verified that the offered strategies 
converge to the AVE (7) under proper selections of the parameters. Theoretical estimation and the 
numerical study of the proposed algorithms indicate that the proposed techniques are applicable for 
determining AVEs. 
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